in 

On 
On 



13 



September 1995 FTUV/95-42 IFIC/95-44 

Invisible Matter 

Lectures presented at XXIII ITEP Winter School, 

Zvenigorod, March, 1995 
(to be published in Surveys of High Energy Physics) 



A.D. Dolgov Q 

Instituto de Fisica Corpuscular - C.S.I. C. 
Departament de Fisica Teorica, Universitat de Valencia 

> 

t-- ! 46100 Burjassot, Valencia, SPAIN 

O 

^ ! Abstract 

in 

These lectures have been given to particle physicists, mostly experimentalists 
Q H - and very briefly and at a pedestrian level review the problems of dark matter. 

6 

H 



The content of the lectures is the following: 1. Introduction. 2. Cosmological 



background. 3. Luminous matter. 4. Primordial nucleosynthesis and the total 



amount of baryons. 5. Gravitating invisible matter. 6. Baryonic crisis. 7. 
Inflationary omega. 8. Intermediate summary. 9. Possible forms of dark matter. 
10. Structure formation: basic assumptions. 11. Structure formations: basics of 
the theory. 12. Evolution of perturbations with different forms of dark matter. 
13. Conclusion. The presentation and conclusion reflects personal view of the 
author that a considerable amount of invisible energy in the universe is in the 
form of vacuum energy (cosmological constant) and possibly in the form of a 
classical field which adjusts vacuum energy to the value permitted and requested 
by astronomical data. 
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1 Introduction 



There are very strong indications that the bulk of matter in the universe is not the normal 
electron-nucleon staff but something unknown which may be made of not yet discovered 
elementary particles, or of some field excitations at astronomically large scales which are 
stable by topological reasons, or maybe even something else. This unknown form of matter 
contributes more than 90% into total mass (or better to say, energy) density in the universe. 
This matter is usually called dark matter. The name reflects the fact that it is not luminous 
and maybe also the level of our knowledge of the subject. A better name is invisible matter 
since it neither emits nor absorbs light. The only observed manifestations of this matter are 
gravitational effects. There are also some theoretical arguments in favor of existence of dark 
matter but, though very persuasive, they are not absolutely compelling. It is the biggest 
challenge in physics of this and may be of the next century to directly observe this new (if 
it is indeed new) form of matter and to study its properties. 
In these lectures I will address he following questions: 

1. Is there indeed dark (invisible) matter and why do we think that it is so? 

2. Could all dark matter be baryonic? If not all, could there still be some invisible 
baryons? 

3. If not baryonic, what is it? There are two ways to address this question: one is to 
use the elementary particle theory for possible candidates for the constituents of dark 
matter and the other is to use the theory of large scale structure formation in the 
universe and the astronomical observations to constraint the properties of dark matter 
particles. 

In the next Section the necessary cosmological background is briefly presented. In sec. 
3 the luminous matter is discussed. Nucleosynthesis constraints on the amount of baryons 
in the universe are presented in Sec. 4. Astronomical data indicating gravitational action of 
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invisible matter are considered in Sec. 5. In sec. 6 the so-called baryonic crisis created by 
a too large an amount of hot gas observed in rich galactic clusters is discussed. Prediction 
of inflationary cosmology on the value of Q is given in sec. 7. In sec. 8 the summary of 
measurements of Q for different forms of matter is presented. A brief discussion of possible 
cosmic relics which may exist in the universe and contribute to the dark matter is given in 
sec. 9. In sec. 10 a critical analysis of the basic assumptions of the theory of large scale 
structure formation is presented. A brief introduction to the basics of the theory of structure 
formation is given in sec. 11. Evolution of density perturbations in the universe dominated 
by different forms of dark matter is described in sec. 12. In Conclusion some speculations 
on the best bet about the form of dark matter is presented and possible experiments and 
theoretical problems, which may shed light on the nature of dark matter, are discussed. 

1 will not review the existing direct experimental searches of dark matter particles in low 
background experiments. It is done by L.Mosca in this School. 

2 Cosmological Background 

The source of gravitational field is the energy- momentum tensor T^ u . In homogeneous cos- 
mological models it is presented in the ideal liquid form which in the rest frame of matter 
can be written as: 

Tfu, = diag(p, -p, -p, -p) (1) 

where p and p are respectively the energy and pressure density. It is usually assumed that at a 
later stage the universe is dominated by nonrelativistic matter when p may be neglected and 
the energy density of nonrelativistic matter, p = p m , is practically equal to the mass density. 
Indeed as a function of the scale factor a(t), which describes the universe expansion, the 
energy density of nonrelativistic matter goes down as p m ~ a -3 , while that of relativistic one 
as p re i ~ a~ 4 . The late dominance of nonrelativistic matter may be not true in models with 
long-lived decaying particles producing relativistic species but in the standard cosmology it is 
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fulfilled with a very good precision. Apart from contributions to T^ v from normal relativistic 
or nonrelativistic matter there may be also a rather mysterious contribution from vacuum 
energy with T^ u = p va c9nu where g^ v is the metric tensor. In contrast to the energy density 
of matter p vac stays constant in the course of the universe evolution. Not long ago a nonzero 
Pvac was considered as rather exotic by majority of astronomers and was usually neglected 
but now with accumulation of new data the attitude is changing. 

Roughly speaking the matter in the universe can be divided into three categories: 1) 
Visible or luminous matter which is observed either by emitted or absorbed light. Almost 
surely this is the normal baryonic staff. 2) Invisible normal baryonic matter. We can make a 
conclusion about its amount considering primordial nucleosynthesis. 3) Invisible nonbaryonic 
matter. This one is observed only by its gravitational action. It can also be divided according 
to the spatial distribution: clustered or uniform. We can "feel" the former by observing the 
motion of the visible tracers like gas around galaxies or galactic satellites while the latter can 
be observed only by cosmological effects and hence the observations are more complicated 
and less secure. In particular vacuum energy belongs to the last case. 

The magnitude of mass or energy density in the universe, p, is usually presented in terms 
of the dimensionless parameter 

n = P / Pc (2) 

where p c is the so called critical or closure density, 

p c = 3H2m Pi = 1.88/4 Q x 1(T 29 g/cm 3 = I0.5hj 00 KeV/cm 3 (3) 
on 

Here mpi = 1.22 x 10 19 GeV is the Planck mass (the Newtonian gravitational constant is 
Gn = m~pf) and h wo is the dimensionless Hubble parameter, 

H = 100/iioo km/sec/Mpc (4) 

The value of H is very poorly known, 0.4 < h wo < 1. The universe age, as we see below, 
favors lower values of H but the recent data indicate that /iioo = 0.7 — 0.8 or maybe even 
higher lU . 



Integrating the Einstein equation: 

\aj om Pl cr 

which governs the expansion of the universe, one can express the universe age through the 
present day value of the Hubble constant and the fractions of energy density of different 
forms of matter (for more detail see e.g. my lecture in the previous ITEP school M or the 
textbooks g |, §, §): 

t * f 1 ^ X (6) 

H JO \/l — Vttot + ^m 1 "' + Qrel x ~ 2 + ^imc^ 2 

where Q m , f2 re z, and f2 mc correspond respectively to the energy density of nonrelativistic 
matter, relativistic matter, and to the vacuum energy density (or, what is the same, to the 
cosmological constant); Q tot = Q m + Q re i + ft vac , and 1/H = 9.8/tJf o x 10 9 yr. 

The universe age is estimated by two different methods: by nuclear chronology which 
uses measurements of the ratios of the long-lived isotopes, l87 Re/ 232 Th or 238 U / 235 U , and by 
the estimated age of old stellar clusters. Both methods give the age in the range of (12 - 20) 
Gyr. The recent analysis gives t u ~ 18 Gyr (for the review and the list of references see e.g. 
paper 0). 

If the universe is dominated by nonrelativistic matter, as is usually assumed, i.e. Q t ot — 
f2 m , then the universe age is approximately given by the expression t u ~ 9.8/ij~ 1 /(l + \/fi Tn /2) 
Gyr. Even with Q = the big universe age and large H are inconsistent. The inconsistency 
becomes stronger for Q tot = 1 predicted by inflationary models. The observed tendency to 
large values of H and t u presents a very strong argument in favor of nonvanishing vacuum 
energy. For Vt tot = Q m + ^vac = 1 eq.(U) gives 

2 1 + V^vac f -s 

t u = t= — In , = (7) 

"&HyVt vac V 1 — lL V ac 

To get, let us say, t u = 14 Gyr for h WQ = 0.75 we need Q vac ~ 0.8. If this is true the 
bulk of energy in the universe is just the energy of empty space, of vacuum. Unfortunately 



our understanding of vacuum energy is very poor. Any reasonable estimate of it gives 
the result which is some 50-100 orders of magnitude higher than the observational limit 
Pvac < 10~ 47 GeV 4 (for the review see refs. ]8|, |Sj). For example there are contributions from 
quark and gluon condensates into p vac which are well established in QCD and which are 
about 10~ A GeV 4 . So there must exist a contribution into vacuum energy from something 
not related to quarks and gluons but with exactly the same magnitude and the opposite sign. 
It is hard to imagine an accidental cancellation with such an accuracy but no dynamical 
mechanism has yet been found. If something is very small, one would naturally expect 
this to be exactly zero and this is a general attitude to vacuum energy. This is one of the 
reasons why models with a nonzero cosmological constant were not seriously considered by 
the establishment. Another reason for that is an unnatural coincidence of the present-day 
values of vacuum energy and the critical energy density of the universe. As we mentioned 
above they are pretty close to each other. However the critical energy density is decreasing 
with the universe age as 1/t 2 while p vac remains constant (at least in the standard model). 
The new astronomical observations changed the attitude to cosmological constant and now it 
is considered more seriously. Moreover, as we see below, models with nonzero vacuum energy 
have some advantages in the description of the structure formation in the universe. To my 
mind the best possibility in solving the cosmological constant problem is the adjustment 
mechanism [[H], [□]] which ensures the cancellation of vacuum energy by an action of a new 
field coupled to gravity or curvature so that its condensate cancels out any initial p vac . This 
cancellation is generically not complete and a noncompensated amount of p vac is always of 
the order of p c . In such cosmological models both vacuum energy and the energy of the 
new field are essential at any stage of the universe evolution (in contrast to the models 
with normal time independent vacuum energy) and this might have an impact on primordial 
nucleosynthesis, structure formation, etc. The progress here is inhibited by an absence of a 
consistent (even toy) model. 
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3 Luminous Matter 



The amount of luminous matter in the universe is estimated in the following way. The flux of 
light, F, from different sources is measured on the Earth and the luminosity of an individual 
source is found by the relation: 

U = AirFtf (8) 

where Zj is the distance to the source. Then one may estimate the average mass density 
of luminous matter as p\ um = J2i Li(Mi/ Li)/V where V is a (large) volume over which 
the averaging is done. Usually the value M/L = 5M /L , which is typical for the stellar 
(galactic) matter, is substituted for the ratio Mi/ Li (here M = 2 x 10 33 g and L & = 
4 x 10 33 erg /sec are respectively the solar mass and luminosity). The measured flux permits 
to find the luminosity density as (see e.g. ref. ||12|| ): 

A = J2 h/V = (2 ± 0.6) x 10 8 h W0 L Q Mpc~ 3 (9) 

i 

where 1 Mpc = 3.1 x 10 24 cm. The result is proportional to the Hubble constant because the 
distance / is determined by the Hubble flow, v = HI and correspondingly l 2 /V ~ /iioo- The 
mass density of luminous matter is given by pi um = X(M/L) = 5(A/L Q )M = (0.7 ± 0.2) x 
10- 31 /iioo g/cm 3 . 

Using these numbers we find for the relative amount of the luminous matter: 

O ta = piujpc = (0.35 ± 0.10)%/Hoo (10) 
This is by far smaller than the amount of nonluminous matter. 

4 Primordial Nucleosynthesis and the Net Amount of 
Baryons 

Though some (or maybe even the bulk) of baryons are invisible there is a way to find their 
total number density with a relatively good precision. This is based on consideration of 
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primordial nucleosynthesis of light elements such as deuterium, helium-3, helium-4, and 
lithium-7. Their calculated abundance are sensitive to the cosmic baryon number density 
or, better to say, to the baryon-to-photon ratio, 7710 = 10 10 ns/n 7 = 10 10 77. The starting 
point of the nucleosynthesis is the fixation of the proton-to-neutron ratio. It is determined 
by the competition between the weak interaction reactions, pe~ <-> nv e and pv e <-> ne + , and 
the universe expansion rate. To the moment when these reactions become ineffective the 
(n/p)-ratio freezes approximately at 1/6. It takes place roughly at T = 0.65 MeV. Later 
on at T pa 0.065 MeV , when the formation of light nuclei begins, this ratio drops down to 
1/7 because of neutron decay. The magnitude of n/p-ratio does not depend upon 77 because 
the corresponding reactions are linear with respect to baryons and so their rate, N/N, is 
independent of the baryon number density. Note that n/p-ratio is rather sensitive to the 
number of different particle species in the primeval plasma and this permits in particular to 
put a bound on the number of different massless (m < 1 MeV) neutrino species. For the 



up-to-date analysis of this bound see refs. |L4|, [15], [U| . 

Since 77 pa 3 x 10~ 10 is a very small number, the amount of the produced nuclei is tiny even 
at the temperatures which are smaller than their binding energy. For example the number 
density of deuterons is determined by chemical equilibrium and is equal to 

where = 2.224 MeV is the deuterium binding energy. One can see that rip becomes 
comparable to n n at the temperature: 

0.065 MeV . , 

Td = : 12 

1-0.03 In r7io 

Above this temperature the number density of deuterons in the cosmic plasma is negligible 
and correspondingly the formation of other nuclei, which proceeds through collisions with 
deuterium is suppressed. This is the so called " deuterium bottleneck" . But as soon as T D is 
reached, the nucleosynthesis goes very quickly and practically all the neutrons, which existed 



in the cosmic plasma at that time, are captured into helium-4. The latter has the largest 
binding energy, Ba He = 28.3 MeV and so in equilibrium its abundance should be the largest. 
Its mass fraction, Yi^He), is determined predominantly by the (n/p)-ratio and is equal to 
2(n/p)/[l + (n/p)} « 25%. 

The time moment corresponding to T = Tjj is determined by the well known relation: 

t D = 0.16 fejpj V2 (^) = 310 sec (^jpj V2 (1 - 0.03 In Vl0 f (13) 

This relation is obtained by equating the expression for the critical energy density of rela- 
tivistic matter 

p c = 3m^/32vrt 2 (14) 

(see eq.(|3|) with H = l/2t) and the energy density of equilibrium relativistic plasma with 
temperature T 

p T = tt 2 ^T 4 /30 (15) 

Here g* is the number of relativistic species in the plasma. For a photon and three types 
of neutrinos with equal temperature g* = 7.25. After electron-positron annihilation the 
temperature of photons became 1.4 larger than that of neutrinos and g* = 3.37. One sees 
that with a larger t] the time to is getting smaller and the number of not-yet-decayed neutrons 
is larger, n n ~ exp(— to/T n ), where r n = 890 sec. Thus Y( A He) is a rising function of r] and 
its variation is relatively weak, with 7710 changing from 1 to 10 the mass fraction of helium 
changes by approximately 10%. Still since Y( 4 He) is rather accurately known it can be used 
as a serious indicator of r] and correspondingly of the total baryonic mass density. 

In contrast to 4 He the fraction of deuterium is a strong and a decreasing function of 
77. Physically it is easy to understand: more baryons are in the plasma, more efficient are 
the processes of disappearance of deuterium through formation of heavier nuclei with larger 
binding energies, like e.g. D + p — > 3 He + 7 and 3 He + n He + 7. The destruction of 
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deuterium in the first process is governed by the equation: 

h D = —avn D n p (16) 

where a is the cross-section of the reaction and v is the velocity of the colliding particles. 
We assumed that the temperature is low enough so that the inverse process is Boltsmann 
suppressed and is not effective. 

Integrating this equation under assumption that r) p = n p / n 7 « 77 = const we find: 

n D = n Di exp[-0. 06r?(3.77 / g*) 1/2 avmpiTi\ (17) 

where nm and Tj are the initial values of the number density and temperature which should 
be taken at the moment when the production of deuterium by photo-dissociation of heavier 
nuclei effectively stops. Eq. ( |17|) illustrates the made above statement about the dependence 
of deuterium abundance on 77. 3 He has the similar dependence on 77 but somewhat less 
strong. The net output of these elements is much smaller than that of 4 ife, D/H ~ 10~ 4 
and 3 He/H (1 — 2) X 10~ 5 by number. 

The fraction of 7 Li depends on 77 nonmonotonically having a minimum at 7/10 ~ 3. 
This more complicated dependence is connected with a competition of different production 
channels. 

Though the abundance of D, 3 He, and 7 Li are not so well known as that of 4 ife, their 
strong dependence on 77 also permits to use them as indicators of the total baryonic mass 
density in the universe. 

Of course the calculations of the light element abundance are not done in this naive way. 
We use it here only for the qualitative description of the phenomenon. The calculations 
are based on the numerical integration of the complete set of kinetic equations with the 
experimentally measured nuclear cross-sections. With the fixed values of the parameters the 
theoretical results are very accurate but the direct comparison with observation is difficult not 
only because of the observational uncertainties but also because of complicated evolutionary 
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effects. Due to the latter the light element abundance in the present-day universe may be 
quite different from the primordial ones. It is relatively simple with 4 He because it can only 
be produced but not destroyed in the course of the universe evolution. The production of 
helium-4 in stars is accompanied by the production of oxygen, nitrogen, and carbon (they 
are called metals by astronomers). Because of that the observations of primordial helium 
are done in the regions of the sky least contaminated by these heavier elements. Moreover 
the data on helium-4 and the heavy elements are linearly extrapolated to the so called 
zero-met allicity state. In this way a rather safe upper bound on Y( A He) can be obtained. 

The evolutionary behavior of other light elements is more complicated, they can be 
both destroyed and produced in stellar processes and the comparison of the theory with 
observations is difficult. Still only two years ago the agreement between the the theory and 
observations was believed to be excellent. In a sense this is true because the theory predicts 
the abundance of different light elements which span the region between 25% for A He down 
to 10~ 10 for 7 Li and this is definitely confirmed by the data. The problem appeared at the 
level of more refined conclusions. The old data for all elements: 4 He, D, 3 He, and 7 Li, were 
in agreement with the theory for the common value of r] 10 = 3.5 ± 1. The best fit to the data 
was obtained with the three known massless neutrinos which was also an argument in favor 
of the overall agreement. With this value of r? the baryonic mass fraction in the universe is 

n B = Ax KTVo/hoo = (1-4 ± 0.4)%^ 2 (18) 
In this case flB/Qium = (2 — 7)hi 00 and the baryons in the universe are mostly invisible. 



A new measurement of deuterium abundance |L7|, [Ll| at high redshift z = 2.9 showed 
a surprisingly high value, D/H = (1.9 — 2.5) x 1CT 4 . It was done by the measurement of 
absorption of light from a distant quasar in a Lyman alpha cloud. The deuterium in this 
cloud did not suffer from evolution and so its abundance should be equal to the primordial 
one. This value is an order of magnitude bigger than the value found in interstellar medium. 
This would request r\ approximately three times smaller than the previously accepted and, 
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if so, practically all baryons in the universe are visible. The analysis of the new data on 
helium-4 also suggests a smaller rj. However another group reported much smaller 
value, D/H = (1 — 2) x 10~ 5 . So the problem remains unsettled. A smaller value of rj creates 
problems with the abundance of 3 He which in this case should be considerably larger than 
that observed in our neighbourhood. For the analysis of the new data and evolutionary effects 



one can address the recent papers [y| [L4|, ^C§ . In accordance with the analysis of helium-4 
made in ref. |TI| the best value for the number of neutrino species is N v = 2.2 ± 0.27 ± 0.42 
if the old value rjio = 3.6 is assumed. Though the central value of N v is below three the error 
bars are generous enough to be consistent with the standard three neutrinos. 

So we have two interesting possibilities which can be resolved with an improved accuracy. 
The first is that rj is close to the old value so that there are invisible baryons and N v < 3. It 
may mean in particular that v T is heavy and unstable. Another more conservative possibility 
would be realized if D/H is high, t] is small and thus all baryons are visible, and N u = 3. 
Small r) may imply a problem with the observed 3 He but evolutionary effects could be 
significant (see ref. fl5fl). A large universe age looks favorable from the point of view 
of primordial deuterium destruction. This can explain a large amount of primordial D/H 
observed in Lyman alpha clouds and a small amount of it in interstellar medium nearby. The 
large tjj in turn is possible only with nonzero vacuum energy. If an adjustment mechanism 
for its cancellation is operating then the noncompensated amount was essential during all 
stages of the universe evolution and this might have an unknown impact on the primordial 
nucleosynthesis. However if we forget about this rather uncertain possibility, the conclusion 
that the primordial nucleosynthesis does not permit the mass fraction of baryons to be larger 
than a few per cent seems to be sound enough. 
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5 Gravitating Invisible Matter 

The observation that there is more gravitating matter in the universe than is directly seen was 



made more than half a century ago by Oort [^TJ for the Galaxy and by Zwicky |22fl for Coma 



cluster. Virial considerations showed that the visible matter alone cannot account for the 
observed velocities. To the present day a very rich data were accumulated which strongly 
support the conjecture that there is a large amount of invisible matter in the universe. 
Velocities of HI gas clouds around galaxies are measured by 21 cm line for hundreds of 
galaxies up to distances r max = 30 Kpc. All this curves becomes flat, v — > const, for r 
outside the luminous centre, r > r ga i = 10 Kpc. The equilibrium velocities are determined 
by the relation 

V - = ^1 (19) 

where Gn = rripf is the Newtonian gravitational constant and M(r) is the mass inside radius 
r. This equation implies that for the mass confined inside r ga i and for r > r ga i the velocity 
goes down with the radius as l/x/r, while for v tending to a constant, M(r) ~ r. Since r max 
is considerably larger than the luminous galactic radius the amount of invisible matter in 
galaxies should be much larger that that of the visible one. 

The law M(r) ~ r implies that the mass density p(r) falls down as 1/r 2 . This takes place 
for the isothermal self gravitating gas as follows from the equations: 

A0 = A-kGnp (20) 

and 

p = p o exp(-0/T) (21) 

There is a very important question how far, that is to what maximum value of r, the 
law M(r) ~ r or v = const is valid. The measurement of HI gas velocities which can be 
done up to 30 Kpc does not show any cutoff and for large spirals give M DM / Mi um = 3 — 5. 
Measurements of velocities of satellite galaxies around the Milky Way and other large spirals 
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can be extended up to 200 Kpc and give Mdm / 'Mi um = 12 — 15. This gives fl sp irais ~ 0.09/iJ" 00 . 
Whether dark halos go beyond this distance is unknown. No cut-off is observed with all 
existing observations. Elliptical galaxies where dark matter is traced up to 100 Kpc by hot 
X-ray gas give M DM /Mi um > 10. An last but not the least (though less secure) the clusters 
of galaxies give M DM /Mi um = 50 — 100. 

The determinations of the gravitating masses by rotational velocities, considered above, 
are valid for gravitationally bound systems and applicable to masses at relatively small scales, 
at most of cluster of galaxies. There is another way to determine gravitating masses which can 
be applied to considerably larger scales. It is based on the measurements of the large scale 
flows. The latter are the so called peculiar velocities i.e. the velocities of different galaxies 
with respect to the general Hubble expansion. By assumption these velocities are induced 
by the density contrast, 6 = (p — p ) / p , where p is the homogeneous cosmological energy 
density. In other words peculiar velocities are proportional to the gravitational acceleration 
created by the contrast of the gravitational potential, 6<p. The theory of this phenomenon is 
simplified because at large scales 6 <C 1 and the linear approximation is valid. 

From the continuity equation, d t p + d(pv) = 0, one can find in the linear regime: 

W = -d t 8 (22) 

where a(t) is the cosmological scale factor. 

For sufficiently large wave lengths (larger than the Jeans wave length, see below sec. 11) 



the behaviour of 6 at matter dominated stage is governed by the equation (see e.g. [^3|, |24j]): 



5 + 2(a/a)5 - ^G N Vtp Q 5 = (23) 

This equation can be easily solved in two limiting cases: 1) Q = 1 for which the rising 
(gravitationally unstable) mode goes as 5 ~ a(t) ~ t 2//3 and 2) Q = for which S = const. 
In the intermediate case, < Q < 1, the solution can be approximately presented as |23|j: 



(6/6) = tt°- 6 {a/a) = HQ 0G (24) 
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So from eq. (|22f) follows 



,0.6 / 6 P\ _ ^n0.6 I 6 P 



W = -HW b P- = --fi U b ^ (25) 

where the biasing factor b is introduced to relate the visible density contrast to the total one. 
Essential assumptions which are made in the derivation of this equation are the linearity i.e 
$tot — fivis/b and the scale independence of this relation. It is also assumed that the visible 
density contrast is given by the excess in galaxy number count 6n/n. All these assumptions 
may be incorrect to a larger or smaller extent. 

It is noteworthy that this equation does not contain the ambiguity connected with the 
cosmological distance scale or (what is the same) with the value of the Hubble constant H. 
Indeed the integration of the equation gives J Hdr ~ Vh where vh is the velocity of the 
Hubble expansion and is directly measurable. Note also that our derivation of eq.(^) is 
oversimplified and can be used for illustrative purposes only; for a more rigorous approach 
see e.g. refs.0, |23| . 

The velocities of galaxies are measured by the red-shift and so only the component in 
the direction to the observer can be found. This extra problem can be solved if the velocity 



field is rotationless, V x v = [25]. It is true for pure gravitational interactions for which 
an initially curlless velocity field remains curlless. Thus the velocity field can be expressed 
through a single scalar function, velocity potential, and this permits to determine Vw from 
observations and to find Q from eq. (pB"D. The recent analysis [^6| of peculiar velocities of 
about 3000 galaxies gives 

n os /b = 0.74 ±0.13 (26) 

So if ft = 1 then b = 1.35 ± 0.23 or if b = 1 then Q = 0.61 ± 0.18. If Q < 0.3 one has to 
introduce antibiasing, b < 1. Though it is not absolutely excluded, the physics of that is 
rather unclear. For the review of this method and possible errors one can address to the 



paper |2"T |. It is argued that the lower bound on independent of biasing is > 0.3. This 



conclusion however is based on the assumption of vanishing cosmological constant. It is 
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worthwhile to redo the analysis without this assumption. Some other potential reasons that 
may invalidate the above conclusion are a possible dependence of b on scale and a relatively 
small size of the analyzed sample of galaxies. Another reason which could also lead to the 
violation of the statement about large Q is a possible nongravitational contribution to the 
large scale flow. If for example the large observed dipole anisotropy of the cosmic microwave 
background (cmb) is not induced by the gravitational acceleration but has its own (intrinsic) 
origin like isotemperature fluctuations with very large (superhorizon) wave length, the results 
based on the measurement of the large scale flow may be invalid. Typically this kind of 
mechanism would create the quadrupole q and higher multipole anisotropy of cmb of the 
same order as the dipole (d) while observationally d ~ 2 x 10 -3 and q = (afew) x 10~ 6 . This 
problem is analyzed in the recent paper []28|| where it is argued that it is still possible to get a 



large dipole while all other multipoles would be small. Note that the model of baryon island 
, where the baryon asymmetry was generated only inside a spherical bubble 



universe 



with the present-day size corresponding to z = 5 — 10, provides q ~ d 2 and very small other 
multipoles. The dipole and other asymmetries in this model are connected with the shift in 
our position with respect to the geometrical center of the island. The smallness of the dipole 
means that we live almost in the Center of the Universe. This is not very natural but rather 
appealing with regard of our vanity. 

6 Baryonic Crisis 



Recent observations f30| , |31|j showed a surprisingly large fraction of baryons in rich clusters 
of galaxies. In the first paper ||30|| , which drew attention to the problem, the mass of the 
hot gas in the Coma cluster was estimated by the intensity of X-ray emission measured by 
Rossat satellite, M gas = (5.45 ±1) x 10 13 /i 10 q M©. It dominates the mass of matter contained 
in the stars, M star = (1 ± 0.2) x l0 13 h{Q M & . The total mass of the cluster was found by 
two methods: by the virial consideration and by the condition of the thermal hydrostatic 
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equilibrium in the gas; both ways give the close results: M total = (6.7 ± 1) x lO l4 hw M Q . 
Thus the ratio of the mass contained in baryons to the total mass of the cluster is 

r^^ = (0.09±0.05)^of (27) 

lvl tot 



Similar results have been obtained in ref. 01 from the analysis of 5 clusters (Coma plus four 
others): 

tb > = (0.03 - 0.08)h w f (28) 



As has been argued in ref. [|32| these results give safe lower bounds on rg and the real 
numbers should be noticeably higher. 

Now if we assume that the mass contained in the clusters represents a fair sample of the 
total mass in the universe then these data together with the nucleosynthesis constraint on 
Qb put the strong upper bound on the energy density in the universe: 

^duster < 0.05r] 10 h^ 2 (0.1/r B ) (29) 

This bound is evidently valid for the clustered matter and is not applicable for the uniformly 
distributed one. 

The authors of ref. |30] have discussed several different explanations of the phenomenon 



like nonstandard nucleosynthesis, nongravitational processes in the structure formation, pos- 
sible mistakes in the standard technique in the estimation of the cluster mass, and an open 
universe model with p tot <C p c . The last possibility is disfavored by the standard inflationary 
scenario (see below). To my mind these data give a good indication to nonzero vacuum 
energy, Q vac ps 0.8. It is uniformly distributed and for the typical cluster size about 1 Mpc 
does not contribute much inside this radius. The nonzero value of vacuum energy is rather 
natural by the reasons mentioned above. Another evident possibility is a large contribution 
from the so called hot dark matter which may be also nonclustered at these scales. This 
however seems to be excluded because the total amount of hot dark matter in the universe 
cannot exceed 30% at least in the models with scale-free spectrum of initial perturbations 
(see below Sec. 12) 
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7 Inflationary Omega 

Aesthetically Q — 1 is definitely the most attractive possibility. It is the only value of Q 
which is not changing with time in the course of the universe evolution. In the universe 
dominated by any normal matter all other values of Q run away from 1. To see this let us 
rewrite eq.(Bl), using definitions (0) and (Bf), in the following way: 



C 



n~ L = i — - (30) 

pa 2 

where C = 3km 2 Pl /87t = const. For nonrelativistic matter with negligible pressure p ~ 1/a 3 
and for relativistic one with the pressure p = p/3, p ~ 1/a 4 . In the general case the behaviour 
of p is governed by the covariant energy conservation: 

p = -3H(p + p) (31) 

and for normal matter with positive pressure p < 1/a 3 when a — > oo. Thus for an open 
universe [C < 0) dominated by nonrelativistic matter Q ~ 1/a — > when a — > oo and for 
the closed universe {C > 0) Q — > oo. 

For vacuum energy p = — p and so p = and pa 2 — > oo. Correspondingly Q — > 1 with 
rising a. It would be very interesting if an adjustment mechanism is found which would lead 
to the equation of state such that pa 2 = const which in turn would lead to Q unchanging in 
the course of expansion. 

Assuming that at the present day Q is somewhere between 0.1 and 2, we see that the law 
(p0|) implies (Q — 1) = O(10~ 15 ) at the nucleosynthesis epoch. At earlier stages the finetuning 
should be even stronger. This is the well known flatness problem which was beautifully solved 
by the inflationary proposal |33|, Typically in most inflationary models the universe is 



dominated by the vacuum-like energy so that p = const and a ~ exp(Hjt). Correspondingly 
(SKnfi — 1) ~ exp(— 2H i t)(Q^ 1 — 1) where r is the duration of inflationary stage, Hj is the 
Hubble parameter during inflation, Q in fi is the value of Q at the end of inflation, and is 



18 



the initial (preinflationary) value of Q. When inflation is over, Q runs away from 1 and at 
the present day reaches the value: 

n 1 - 1 = (nj£, - i)( Zinfl + iy 2 ( Pmfl / Po ) » (nrj, - i)(t r /t ) 2 (32) 

where subzero means the present-day values, (z in fi + 1) = (ao/o^/j) ~ (Tr/T ) is the red- 
shift corresponding to the end of inflationary era, Tr is the temperature of the universe 
heating at the end of inflation, and To = 2.7 K is the present-day temperature of the cosmic 
microwave background radiation. We assumed that the energy density at the end of inflation 
is Pinfi = (tt 2 /30)g^AT A where gf is the number of relativistic species at T = Tr and A > 1 
is the coefficient which accounts for the slow process of heating due to the weak inflaton 
coupling. The last equation in (|32| ) was obtained under assumption Ag* (p C mb)/po ~ 1- The 
necessary duration of inflationary stage is approximately equal to Hit ~ ln(Tft/To) = 52 + 
ln(TR/10 10 GeV). This gives CIq 1 — 1 = 0(1). Even a slightly longer inflation would result in 
| Qq 1 — 1 |<C 1 and a slightly shorter one would be in a gross disagreement with observations. 
Typically inflationary models give Hjt 3> 1 which results in | f^ 1 ~ 1 1^ 1 i 10~ 4 . The 
deviation from 1 in the r.h.s. of this relation is connected to the density inhomogeneity 
at the present horizon scale. So Q = 1 is considered as a robust prediction of inflationary 
cosmology. 

Recently with accumulation of new astronomical data which may request Q < 1 there 
appeared models where a possibility of having Q ^ 1 in inflationary frameworks was advo- 
cated P5]| . All these models request a tuning of inflationary and postinflationary stages in 



such a way that exp(2#r) = (Tr/T ) 2 (^o ^l) -1 . This condition looks very strange because 
physically these epochs are not related. The condition of small density perturbations which 
is not easy to realize in these models imposes an extra restriction making these models even 
less natural. 

Thus even if Q ^ 1 is in principle possible in inflationary models the flatness problem 
in such models is similar to that in the old Friedman cosmology. Still it is very difficult 
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(if possible) to create our universe without inflationary stage and if it is indeed proven that 
^matter < 1 a much better solution which naturally agrees with inflation is a nonzero vacuum 
energy such that Qtot = ^matter + &>vac = 1- On the other hand if we observe that VL to t ^ 1 it 
would be a strong argument against inflation (1 have to admit that this opinion is not shared 
by everybody). 

8 Intermediate Summary 

Let us summarized what is known about energy density of different forms of matter in the 
universe. 

1. The luminous matter which is undoubtedly made of the usual particles, protons, 
neutrons, and electrons, contributes very little to the total mass of the universe: 
^lum ~ 3.5 X 10 3 /iioo- 

2. Primordial nucleosynthesis gives for the total amount of baryons = 4x 10 -3 77io/i7oo 
with r/io = 1 — 7. At the moment there is no consensus whether 77 is relatively high 
(near the upper bound) or at the lower end of the permitted interval. In the first case 
we should expect plenty of invisible baryons in the universe while in the second all 
baryons may be visible. 

3. Measurements of rotational velocities in gravitationally bound systems give, depending 
on scale, Q = 0.1 — 0.3. 

4. From the large scale flows follows Q; s / = [6(0.74 ± 0.13)] 5 / 3 with the biasing parameter 
6 = 0(1). 

5. From the lower bound on the universe age tu > 12 Gyr and h wo > 0.7 follows Vt matter < 
0.2. 
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— 1 /2 

6. Measurements of X-ray emission from rich galactic clusters give ^clustered < 0.15/i 10 o 
0.55h% 2 ). 

7. Inflationary cosmology predicts fltot = 1 with a very good accuracy. 

It is rather difficult to say now what is the best choice in this situation. The answer 
would be different depending on the personal prejudices of the author. To my mind the odds 
are the best for the following: flf, ar « fl vis « 0.5%; Qdm = 0.1 — 0.2, and Q vac = 0.8 — 0.9 
with Vt tot = 1. 

9 Possible Forms of Dark Matter 

The simplest and most economical is to assume that all dark matter consists of the usual 
known particles: protons, neutrons, and electrons, without unnecessary introducing quan- 
tities according to Occam. Unfortunately this idea of a large amount of dark baryons con- 
tradicts the nucleosynthesis bound and the only visible way to overcome this bound is to 
conceal the majority of baryons in primordial black holes. It can be done but with a rather 
exotic model of baryogenesis which gives a very big baryon asymmetry inside relatively small 
bubbles [ 36] . When the temperature drops below the proton mass, very big density perturba- 



tions evolve which result in almost 100% collaps of the regions with high baryon asymmetry. 
As is argued in ref . ||36|| (see also if baryonic black holes have been formed sufficiently 
early the bound ( f[8~D may not be applicable and Qb close to 1 is permitted. 

There are many direct astronomical restrictions on the baryonic dark matter in all possible 
forms which have been reviewed in lectures |38|] . If the dark baryons are in the form of black 



holes the following limits are to be fulfilled. Black holes with M > 10 6 M o would lead to 
formation of too heavy galactic nuclei. Black holes with M > 10 4 M Q would disrupt globular 
clusters. Black holes with M < 100M Q , which have been formed in the process of heavy 
star evolution, are excluded because the star-progenitors would eject too much metals into 
interstellar medium. However primordial black holes are not subject to the last bound and for 
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them M < 1OOM is permitted. Moreover, as we mentioned above, even the nucleosynthesis 
bound for them may be avoided. 

Another possible candidate for dark baryons are brown dwarfs, light star-like bodies 
consisting of hydrogen and helium with the mass M < O.O8M , so that the reaction of 
hydrogen fusion is not ignited. Dark baryonic matter could also be in the form of the so 
called "Jupiters", planetary-type objects with M < 1O _3 M . Compact stellar remnants: 
white dwarfs and neutron stars (M = (0.4 — 2)M ) are also viable candidates. And at last 
one can expect some dark baryons to be in the form of the cold diffuse gas clouds. 

During the last 3 years several groups of astronomers (EROS, MACHOS, OGLE) are 
looking in our galaxy for dark star- or planetary-size objects with the mass in the range 
(10~ 8 — 1O 3 )M using the gravitational microlensing effect f3~9fl . These objects have got the 
name MACHOS which is the abbreviation for Massive Astrophysical Halo Objects. When 
such an object happened to be on the line-of-sight between the observer and a star the 
brightness of the star may considerably increase due to gravitational focusing of light. The 
luminosity curve should have a symmetric form as a function of time and should not depend 
upon the light frequency. Around 100 events of this type have been observed to the present 
time. Unfortunately it is difficult to judge if these events are induced just by the ordinary 
stars or by the invisible celestial bodies which are looked for. The experiments are now in 
progress and more data are being accumulated. Moreover the new series of observations are 
now done in real time when the change of star brightness is immediately registered, while 
previously the analysis of the accumulated data was performed. So hopefully in the nearest 
future a more decisive conclusion can be made. Still even if all the observed events are 
induced by MACHOS is seems very plausible that they cannot account for all dark matter 
in the universe. For the recent review see e.g. ref.|40| . 

If we confine ourselves to the particles which are known to exist, the next simplest 
possibility for dark matter are massive neutrinos with m = O(10 eV). Though we cannot 
predict theoretically the magnitude of the neutrino mass, there is no known principle that 
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forbids nonzero m v and so it is quite natural to expect them to be massive. Unfortunately 
the theory of large scale structure formation excluded massive neutrinos as dominant bearers 
of dark matter at least in simple models with scale free spectrum of initial perturbations. 
It is an interesting question whether there exists a spectrum of initial perturbations which 
permits to describe the universe structure with neutrinos and baryons only. 

If baryons and massive neutrinos cannot account for the observed Qdm there are still 
plenty of possibilities but now we have to turn to hypothetical not yet discovered particles 
or new objects or fields. The most popular now candidate for the dark matter particle is 
the lightest supersymmetric particle (LSP). In simple supersymmetric models this particle 
is stable due to conservation of R-parity, R = (— \^ B+L+2S where B and L are respectively 
baryonic and leptonic charge and S is the particle spin. So defined -R-parity has the value 
(+1) for all usual particles and (—1) for all (not yet discovered) superpartners. 

In models with low energy supersymmetry the SUSY breaking scale is phenomenologically 
acceptable in the interval Asusy — (0-1 — 1) TeV. With this scale of supersymmetry breaking 
the theory naturally predicts the relic abundance of LSP's near Q = 1. Indeed the number 
density of the primordial LSP's which survived to the present epoch can be roughly estimated 
0: 



as 



n L sp = — (33) 

crvm LSP mpi 

where a is the cross-section of annihilation of LSP's, v is their velocity, and n 7 = 400 /cm 3 
:r d 

where a ~ 10~ 2 and m LS p ~ A S usy we get 



is the number density of photons in cosmic microwave background. With av = <y 2 /m 2 LSP 



Q LSP « 0Mh^ (m LSP /TeV) 2 (34) 
More details about and a review on supersymmetric dark matter can be found e.g. in paper 
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However the conservation of R-parity is not obligatory. No theoretical principle like e.g. 
gauge invariance demands that. So it is quite natural to expect that it is indeed broken. 
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Moreover if we believe the principle, which has been established by the development of the 
elementary particle physics during the last 30 years, that everything that can be broken 
should be broken, R-parity is surely broken. If this is the case the LSP should be unstable. 
Still depending on the mechanism and the strength of breaking the life-time of LSP's can 
be very long like e.g. the life-time of proton which is stable in all known experiments and 
observations despite an almost certain baryonic charge nonconservation. If so, the long-lived 
LSP can well be the dark matter particle. If however their life-time is cosmologically short, 
they could not contribute to dark matter and other candidates are wanted. 

In models with broken R-parity an interesting candidate for dark matter particles can be 
massive majorons [[43], (PJ i.e pseudogoldstone bosons appearing due to spontaneous breaking 



of leptonic charge conservation. The majorons may have mass in KeV range and the number 
density corresponding to Q = 0(1). The universe structure formation in such model has 
been considered in ref . [ i~4" and reasonably well agrees with the observed picture (see section 



12). 

There could exist some other stable (or cosmologically long-lived) particles that may be 
bearers of dark matter. They could be produced thermally, like LSP's or (massive) neutrinos 
in hot equilibrium primeval plasma. Depending on the strength of their interactions and the 
magnitude of the mass they may be relativistic when they are decoupled from the cosmic 
plasma (e.g. neutrinos) or nonrelativistic (e.g. LSP's). Another possible mechanism of 
production may be related to a phase transition and in this case even very light particles 
would be produced in nonrelativistic state. An example for that is the well known axion which 
is also considered as a serious candidate for dark matter particle. The effective potential for 
axion field has the form of Mexican hat after the underlying [/(l)-symmetry is broken at 
high temperature. Later on at QCD energy scale the bottom of the Mexican hat potential 
became tilted and because of that axions acquire a nonzero mass. At that moment coherent 
oscillations near the potential minimum begin. This can be understood as a production of 
Bose condensate of axions (at rest). Though the mass of axion is extremely small, m a = 
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O(10~ 5 eV), the axions, thus produced, remain nonrelativistic because of negligibly weak 
interactions with the plasma. 

Some other, maybe more exotic, and because of that more interesting, candidates for 
dark matter are considered in the literature like e.g. topological or nontopological solitons 
with astronomically large size. Topological defects (or topological solitons) might be formed 
in the course of phase transitions accompanying spontaneous symmetry breaking when the 
universe expanded and cooled down. Depending upon the symmetry group different kinds 
of objects can be created. If the corresponding symmetry group is U(l) then vortex lines 
(or cosmic strings) would be formed which could both make contribution into dark matter 
and create density inhomogeneities |36|. If a discrete symmetry is broken domain walls 
separating different vacua are created. Domain walls with a natural (whatever it means) 
scale are cosmologically forbidden since they would create too large angular fluctuations of 
temperature of cosmic microwave background [17] but if the symmetry breaking scale is 



extremely low, they are allowed and can even serve as seeds for structure formation [4£ 



Spontaneous breaking of 577(2) or 0(3) gives rise to production of monopoles. These three 
types of objects exhaust all possible topologically stable solitons in three dimensions. There 
could also be unstable but long-lived ones which might be cosmologically interesting. For 
the review see e.g. refs. [H9|, 



I think that a very interesting possibility for dark matter is a classical field which possibly 
compensates nonzero vacuum energy. The theory of such field has not yet been seriously 
considered (even at a toy model level) but the cosmological constant problem makes the 
existence of such a field very probable. 

At the present level of our knowledge we cannot say what is the theoretically preferable 
candidate for dark mater particles. We only can make a judgement on the basis of the 
universe structure formation comparing the picture predicted by a model with a particular 
kind of dark matter particles and astronomical observations. 
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10 Structure Formation. Basic Assumptions. 

The universe is believed to be homogeneous on the average at very large scales, I > l\ ~ 
100 Mpc. This hypothesis is based on the observed isotropy of the cosmic microwave back- 
ground and on the absence of the observed structures at I > l±. However the present-day 
accuracy of measurements is not good enough to exclude structures at larger scales so the 
last statement may be questioned. At smaller scales the universe is very inhomogeneous: 
there are galaxies, their clusters and superclusters, large size voids, and possibly even a 



periodicity (or at least structures) with the characteristic wave length about 100 Mpc f51 ]. 

The task of the theory of large scale structure formation is to give a qualitative description 
of the distribution of matter in the universe. This is done with the following four essential 
assumptions of a different level of creditability: 

1. It is assumed that the structure of the universe has been formed as a result of evolution 
of initially small fluctuations under the action of gravitational forces. This is the only solid 
assumption of the theory with which practically everybody agrees. Still even this sometimes 
is questioned and speculations about possible nongravitational processes or initially large 
fluctuations are put forward. 

2. The next very important assumption concerns the nature and the spectrum of initial 
perturbations. For many years the origin of primordial fluctuations at very large scales re- 
mained mysterious. A natural idea that density inhomogeneities were generated by quantum 
or thermal fluctuations could not be realized in frameworks of the Friedman cosmology with 
the power law expansion regime, a(t) ~ t a , because the corresponding wave length were 
negligibly small in comparison with the necessary astronomical scales. This difficulty was 
resolved in inflationary scenario where the wave lengths of perturbations were expo- 
nentially inflated, Iq — > Zoexp(ifr), and what's more their amplitudes were amplified in the 
course of expansion. In a sense the plan was overfulfilled because the density perturbations 
generated in this way would be too high unless an unnaturally weak interaction of the infla- 
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ton field is assumed. It can be shown that for the inflaton with the self-interaction potential 
U(4>) = A0 4 the magnitude of density perturbations with wave length I at the moment of 
horizon crossing is roughly given by the expression: 

5 = J. « 100VX[1 + 0.01 \n(l/l gal )f 2 (35) 



For the details see either the original papers [5l| or the book fl5"3f . Note that to create 



the necessary magnitude of fluctuations, 5 < 10 -4 , bounded by the data on the angular 
fluctuation of the background radiation temperature, one needs A < 10 -12 . 

A competing source of density perturbations may be astronomically large topological 
defects mentioned in section 9. The simplest and the most attractive possibility to my mind 
are cosmic strings. Unfortunately the accepted SU(3) x SU(2) x [/(l)-theory of particle 
interactions does not lead to cosmic strings (neither to any other mentioned above stable 
topological defects) and at the present level of our knowledge we cannot say if higher energy 
modifications of the theory definitely request existence of such objects. 

If we neglect the log-dependence on the scale, the spectrum of fluctuations given by 



eq . (p5|) , coincides with that proposed by Harrison |54| and Zeldovich |55| some time before 
inflationary cosmology was discovered. This is the so-called flat or scale-free spectrum of 
fluctuations. It is the simplest possible primordial spectrum which does not introduce any 
particular scale to the theory. With an advent of inflation this kind of spectrum got the- 
oretical justification and now is commonly used as the basic spectrum in calculations of 
structure formation. Sometimes as a simple generalization an arbitrary power law spectrum 
is considered: 

The Harrison-Zeldovich spectrum corresponds to n = 1. The spectrum with n / 1 may 
appear in some inflationary models |56j but with n not much different from unity. If we 



permit in principle an existence of primordial spectrum with an arbitrary n, introducing in 
this way a new scale to the theory, a combination of several terms with different n as well as 
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a more complicated function are also permitted but without a guiding principle for choosing 
a particular form of the spectrum, the theory would completely loose its predictive power. 

3. Because of stochastic nature of the fluctuations one should consider average values 
and this explains in particular the consideration of 5 2 instead of 6 in eq . (|3~6D . Stochastic 
properties of S(x) are fixed by the following hypothesis. It is assumed that the fluctuations 
are gaussian with delta-correlated Fourier modes. In other words if the relative amplitude 
of the fluctuations is Fourier expanded: 

S(x) s m = J l^ffle* (37) 

the averaged product of the mode amplitudes satisfies the condition 

(5(k)5(q)) = (^) 3 f(kmk-q) (38) 

(hopefully the delta-function here is not confused with the relative magnitude of density 
fluctuations S(x) and its Fourier transform 5(h)). 

It is also assumed that the spectrum is isotropic so that the function / is a function of 
the absolute value of k only, f(k) = f(k). (Note in passing that this assumption is not 
obligatory but we do not know any reasonable physical mechanism for the anisotropy.) 

Using these equations we can easily find 

(5 2 (x)} = J d 3 kf(k) (39) 

With f 2 = const /k n+2 it coincides with eq. fl36|) . 

4. Except for the three hypothesis discussed above, the structure formation evidently 
and heavily depends upon the properties of matter which presents building blocks for the 
structures. In contrast to the previous three assumptions there is no consensus with respect 
to the universe chemical content and a plethora of different possibilities is proposed. From 
the point of view of structure formation all possible forms of matter particles are classified 
as cold, hot, and, for the intermediate case, warm. Cold matter or CDM (here "D" stands 

28 



for "dark" because it is assumed that the universe is dominated by dark matter) is formed 
by particles which were nonrelativistic at a rather early stage when the galactic size crossed 
the horizon. Hot matter (or HDM) is formed by particles which were relativistic at that 
moment, and correspondingly warm particles were semirelativistic. This property of dark 
matter particle determines the characteristic size of objects formed. 

At the present day there are several models which more or less successfully describe 
structure formation under the first three assumptions and with different hypotheses about 
the matter content of the universe: CDM plus vacuum energy mixed hot plus cold 
dark matter |58|], unstable particles decaying into relativistic species plus CDM U59| , HDM 
plus astronomically large seeds like topological defects. These models request comparable 
contributions of different forms of matter which are not physically related at least at the 
present level of our understanding and by this reason do not look natural. On the other 
hand we already have an example of such unnaturalness realised in the universe, namely the 
energy density of baryons is relatively close to that of dark matter (compare Qf, ar with £Idm) 
though they may be different by several (many!) orders of magnitude. At the moment it is 
not clear if this rough coincidence is an indication of an unknown deep physics, or there is 
something wrong in our understanding of dark matter, or we have to invoke the anthropic 
principle to explain that. Models with a single dominant form of matter demand a deviation 
from the flat spectrum of primordial density fluctuations. 

It is usually assumed that the dark matter is collisionless (or in other words its selfinterac- 
tion as well as interactions with other particles may be neglected) but this is not necessarily 
so and the models with self-interacting dark matter are also considered [61], |44| . 

We will briefly discuss the relevant properties of models of structure formations with 
different forms of dark matter particles or fields below. 



29 



11 Structure Formation. Basics of the Theory. 

The evolution of small density perturbations is described by the usual hydrodynamic equa- 
tions of selfgravitating perfect fluid with some complications related to the expanding cos- 
mological background and effects of general relativity. The dynamics of liquid element is 
governed by the forces of gravity and pressure. If the pressure dominates then the pertur- 
bations evolve as sound waves. For the collisionless liquid like e.g. that of neutrinos the 
perturbations are erased by free-streaming which formally very much resembles the effect of 
pressure. In the case when gravity dominates, the element of liquid would be indefinitely 
compressing if the equation of state remains unchanged and correspondingly no new pressure 
forces come into play. In this case the density contrast would be rising and this results in 
formation of structures in the universe. Evidently with an increasing size or mass of a system 
gravity would ultimately dominate pressure while for small systems pressure is normally the 
dominant force. The boundary value of the mass and size are called respectively Jeans mass, 
Mj, and Jeans wave length, Xj. They are named so after Jeans who was the first to study 
the phenomenon of gravitational instability at the beginning of this century. 

The basic equations which describe evolution of density perturbations in the nonrela- 
tivistic perfect fluid with mass density p and pressure density p are the continuity equation: 

d tP + 8{pv) = 0, (40) 

the Newtonian equation of motion (mx = F): 

d t v + (vd)v + dp/p + d(f) = 0, (41) 

and the Newtonian equation for the gravitational potential (j>: 

A0 = AnpG N (42) 

Perturbative expansion of these equations over homogeneous background p = po + pi(t,x), 
v = vi (x, t) , etc, permits to derive the second order linear differential equation describing 
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evolution of density perturbations: 

6 - v 2 s A5 - 4np G N 5 = 



(43) 



where v s is the velocity of sound. 

Making Fourier transform of this equation one sees that short wave modes are oscillat- 
ing while the long wave ones are exponentially rising, 5 ~ exp 

(^ttGjvPo - v 2 s k 2 1). The 

boundary value of the wave vector k 2 = AttGnPo/v 2 is called the Jeans wave vector and the 
inverse quantity Aj = 2n/kj is the Jeans wave length. The Jeans mass is the mass contained 
inside the volume bounded by Aj and is equal to Mj = 47r(Aj/2) 3 p /3 = 
In realistic cosmological situation these expressions remain true with the substitution the 
cosmo logical energy density po(t) ~ 1/t 2 instead of po — const . Because of that the charac- 
ter of instability becomes different, it is no more exponential but of a power law as we have 
seen in sec. 5. 

To derive the equations of motion in Newtonian approximation with the account of 
the universe expansion one has to change all the derivatives to the covariant ones in the 
Robert son- Walker metric, to substitute as zero order approximation the cosmological time 
dependent energy density po(t) (as we have already mentioned) and the nonvanishing zero 
order velocity, Vq = Hf. Otherwise the procedure is exactly the same as above. The resulting 
equation for the Fourier transformed amplitude of density perturbations 5k is similar to eq. 
(|43|) with an extra "friction" term proportional to H = a/a: 

8\ + 2(d/a)4 + (v 2 s k 2 /a 2 - AnG N n Po )8 k = (44) 

The limit of small kja of this equation gives eq. (p3|) . More details about the material 
presented above and below can be found in the books KB], O or in the lectures p§, W\ . 



It is essential for structure formation that perturbations in nonrelativistic matter do not 
grow in the universe dominated by a uniformly distributed relativistic matter. In this case the 
term which drives the instabily in eq. fl44]) is changed to 47ip G]ye where e = p m /po, p m is the 
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energy density of nonrelativistic matter, and po is the dominant cosmological energy density 
of relativistic matter. Since by assumption e < 1, the driving force becomes negligible and 
rising modes are not generated. If the relativistic matter is itself inhomogeneous then there 
are unstable modes but usually inhomogeneities in relativistic matter are quickly erased by 
free streaming or diffusion. Correspondingly cosmological inhomogeneities started to grow 
only at matter dominated stage. 

General relativistic treatment of the problem is more complicated. One has to solve the 
Einstein equations for the metric with corrections from inhomogeneous fluctuations. There 
is an ambiguity related to the freedom in choice of different coordinate systems (the so called 
gauge freedom), so that one has to distinguish between the gauge artifacts and real physical 



effects. This analysis was first done by Lifshits in 1946||63||. The gauge independent formalism 



was developed by Bardeen |64| relatively recently. In many cases however it is sufficient to 
use the Newtonian approximation for qualitative (and often quantitative) understanding of 
the phenomena. 

While the analysis of the evolution of small perturbations is relatively simple the non- 
perturbative regime which started when 5 = 0(1) is much more complicated. No reliable 
analytical or semianalitycal methods have yet been developed and the calculations are made 
by numerical simulations with a large number of noninteracting (except for gravity) parti- 
cles. By necessity many essential properties of the system are neglected but nevertheless one 
may hope that this approach gives a reasonable description of the structure formation at the 
stage of large perturbations. 

12 Evolution of Perturbations with Different forms of 
Matter. 

Let us first consider density perturbations in the cosmological gas of light neutrinos with 
mass in 10 eV range. At an early stage when neutrinos are relativistic density contrasts 
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are smoothed down by their free motion. The perturbations could start rising only when 
neutrino momentum is redshifted down to nonrelativistic values. The characteristic size of 
the fluctuations which survived the erasure by free streaming, is given by the path which 
neutrinos could travel till they became nonrelativistic, l m ~ 2t m . Here t m is the universe 
age at that moment. It is assumed that neutrinos moved with the speed of light and the 
coefficient 2 is connected with the universe expansion (in the static universe the path would 
be just I = t). We find t m using equations (|T^JT5| ): t m = 4 x 10 20 cm(10 eV/T m ) 2 . We choose 



T m = because the average neutrino energy in thermal equilibrium state is (E) « 3T. 

To calculate the magnitude of l m at the present day we have to multiply it by the red shift 
z m + 1 = T m /T° where T° = 1.93 .fT is the present day temperature of the cosmic massless 
neutrinos. Thus we get 

= 16 Mpc(10 eV/T m ) = 48 Mpc(10 eV/m v ) (45) 

The mass of matter inside the radius l m /2 is M m « 3 x 1O 15 M (1O eV/m u ) 2 . It is 3-4 orders 
of magnitude bigger than galactic mass. Note that 1$ and M m coincide correspondingly 
with Xj and Mj given above if we formally take the speed of sound equal to the speed of 
light. This is indeed true for relativistic gas. Eq. (|45|) gives the mass of the smallest size 
objects which could be initially formed in neutrino dominated universe. Smaller objects is 
difficult to create and this was the reason why a beautiful idea of neutrino universe was 
abandoned. Particles with a larger mass like e.g. m = 300 eV would create structures with 
characteristic size close to that of galaxies. As we mentioned above these are particles of 
warm dark matter. There is a renewed interest to these particles now. Unfortunately no 
particles with such mass are observed experimentally but there are some theoretical models 
(see e.g. refs.[(|3|, [HJ]) where such particles may exist and contribute to dark matter. 



During several last years, before the COBE measurements |)5] of the large angle anisotropy 
of the cosmic microwave background radiation, the model with cold dark matter was def- 
initely the most popular (for the review see e.g. refs. |58|, |66|). Based on a very simple 
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assumption of scale- free spectrum of primordial fluctuations (Harrison-Zeldovich spectrum), 
the model reasonably well described the galaxy distribution at the scale of tens megaparsec 
with the only free parameter, the overall normalization of the spectrum. However the COBE 
measurements have fixed the normalization at large scale end of the spectrum and with this 
normalization the prediction for smaller scales became a factor two above the observations. 
An evident possible cure is to change the shape of the spectrum of primordial fluctuations. 
It seems to be a very interesting possibility but since this idea is not supported by the in- 
flationary scenario, the major line of investigation is the consideration of different forms of 
dark matter with the same flat spectrum of density perturbations or maybe an introduction 
of the cosmological constant. The basic idea of all these models was to suppress the power 
of the evolved spectrum at smaller scales relative to that at larger (COBE) scales. 

One possibility is to shift the epoch of matter dominance (MD) to a later stage in a simple 
cold dark matter model. Since the characteristic scale at which perturbations started to rise 
in this model is determined by the horizon size at the onset of the MD stage, shifting it to a 
later moment gives less time for rising of the fluctuations and correspondingly less power at 
galactic and cluster scales. This goal can be achieved if one assumed that universe is open 
so that hf 00 Q pa 0.2. This model is disfavoured by inflationary scenarios which (at least in 
simple versions) predicts Q = 1. One can recover this prediction of inflation in the universe 
with low matter density if the cosmological constant A (or in other words, vacuum energy) 
is nonzero. As we have mentioned in sec. 2, the recent data indicating a rather high value of 
h 1Q0 pa 0.8 support the idea of nonzero A with the fraction of the vacuum energy Q vac = 0.8. 



The models with nonzero A give a satisfactory description of the observed structure |57j with 
flat spectrum and COBE normalization. 

A mixed (hot+cold) dark matter scenario can also do the necessary job of diminishing 
the power at small scales because (initially flat) perturbations in hot dark matter are erased 
at scales smaller than ~ 1O 14 M by free streaming if the dark matter is collisionless as is 
the case of neutrinos. A good description of the structure requests 70% of CDM and 30% of 



34 



HDM [58]. It gives even better description if there are two equal mass neutrino species each 



with the mass 2.5 eV |67| as is suggested by the recent indications of neutrino oscillations 
by the Los Alamos group. 

Recently there appeared a renewed interest to the idea of structure formation with unsta- 
ble particles It is assumed in these models that there exists a massive long-lived particle, 
usually tau-neutrino with the mass in MeV range which decayed into massless species at the 
epoch when the mass density of the parent particles dominated the energy density of the 
universe. Correspondingly the present-day energy density of relativistic particles would be 
bigger than in the standard scenario and the onset of MD stage would take place later. 

The common shortcoming of these models is that they all demand a certain amount 
of fine-tuning. Generally one would not expect that the contribution from hot and cold 
dark matter into the universe mass are about the same, they may differ by many orders of 
magnitude. One would also suspect that the vacuum energy which remains constant in the 
course of the universe expansion is by no means related to the critical energy now (which 
goes down with time as m 2 Pl /t 2 ) Q. The models with unstable particles mentioned above 
are also based on the assumption of two independent components: the massive unstable 
particles themselves and unrelated cold dark matter. This is definitely unnatural and this 
shortcoming stimulated search for other models. Recently in ref. a return to to the 
model with a single cold dark matter component universe (of course except for baryons) was 
advocated. For successful description of the observed structure the authors need the power 
index of the spectrum n = 0.8 — 0.9, a low value of the Hubble constant, h wo = 0.45 — 0.5 (to 
be compatible with large Q), and a large contribution of tensor perturbations (gravitational 
waves) into quadrupole fluctuation of the background radiation temperature measured by 
COBE, C T /C S = 0.7. 

Another attempt to overcome unnaturalness of multicomponent dark matter model has 

2 Adjustment mechanism [|o[ |ll| (for more references see reviews || ||] and the recent paper |||) however 
generically gives the noncompensated amount of vacuum energy of exactly necessary order m 2 Pl /t 2 and so 
the models with vacuum energy may be more " equal than others" . 
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been done in paper ||44j| . In this paper a model is considered in which unstable particles and 
the particles of cold (or possibly warm) dark matter are closely connected. In fact the decay 
of the former produces particles of the present-day dark matter. A necessary background 



model of this kind in particle physics was proposed some time ago [|7C], |7T[ as an attempt 
to find a phenomenologically acceptable description of R-parity breaking. The underlying 
mechanism is the spontaneous breaking of leptonic charge conservation at electroweak scale. 
The model contains a Majorana type tau-neutrino with the mass around MeV which decays 
into massive but light Majoron J with mass in KeV region. Life-time with respect to this 
decay, as estimated in refs. |70], could be of order days or years depending on the values of 



parameters, i.e. just in the interesting for us interval. It is worth noting that there is no stable 
SUSY particle in this model so the dark matter cannot be associated with it but the model 
itself produce a candidate for dark matter, namely a massive Majoron. The model possesses 
some features like sufficiently large diagonal coupling of v T to majorons or selfinteraction of 
majorons which rather naturally permit to resolve appearing cosmological problems in par- 
ticular the problem of extra massless particle species during primordial nucleosynthesis. The 
cosmological properties of this model are rather unusual and are interesting by themselves. 
Dark matter particles (majorons) in this model are strongly self-interacting and thus the 
structure formation in this model is different from the traditional one with collisionless dark 
matter particles |X| |6l], [72|, |44|]. In particular the shape of galactic halo would depend on the 



dark matter selfinteraction. The recent data [73| might give an indication that collisionless 
dark matter gives a poor description of the shape of halo in dwarf galaxies. 

13 Conclusion. 

It is difficult to write a conclusion on something which is so evasive as dark matter. Though 
it is practically certain that dark matter exists, it is unknown what it is and even if there 
are one or more forms of dark matter. So at the present level of our knowledge we can only 
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make bets which are very little constrained by the poor available data. 

One of the most important questions is whether fl is indeed equal to 1 as predicted by 
inflation or it is considerably smaller. In the last case the universe age problem is not so 
severe with a large Hubble constant but still if foioo > 0.75 and tu > lAGyr even Q = does 
not help much. It is very important to know H more accurately since its value may give a 
weighty argument in favor of nonzero cosmological constant A. With nonzero cosmological 
constant Q ma tter could be relatively small (while Vt tot = Q ma tter + ^vac = 1) and the universe 
age problem does not exist. Structure formation in models with A ^ looks reasonable 
too and the unnaturally close values of p c and p vac may be understood with an (unknown) 
adjustment mechanism. So one of the most important problems is to find a workable model 
for the latter or in more general terms to solve the cosmological constant mystery. If it is 
indeed solved by an adjustment mechanism one should expect that Sp vac ~ p c during all 
the history of the universe. Hence the cosmology with adjustment of vacuum energy may 
give somewhat different predictions for light element abundance, structure formation, etc. 
Probably the baryonic crisis is also resolved by a nonzero cosmological constant. 

Speaking about more traditional forms of dark matter we have to note that nonzero p vac 
is not sufficient to solve all problems with dark matter and some amount of "normal" dark 
matter is necessary. What is it is a big question. It may be LSP's or axions, or massive 
neutrinos. Though the latter are disfavored by the theory of structure formation, it may be 
premature to exclude them completely as a dominant dark matter (together with vacuum 
energy). If R-parity is not conserved LSP's are probably excluded but a massive majoron 
in this case becomes a very interesting candidate. Accelerator experiments in search of 
supersymmetry are very important for resolution of this problem. 

Invisible baryons may also exist and to this end a clarification of situation with primordial 
nucleosynthesis is very desirable especially a resolution of the ambiguity with the observation 
of primordial deuterium. In particular it may settle down the problem with the number of 
massless particle species present at nucleosynthesis and exclude or confirm a hypothesis of 
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MeV-tau-neutrino. 

In a whole we can expect a really exciting development in the nearest future. 
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